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Abstract 

X 

ti- 
lt is known that the Kadomtsev-Petviashvili (KP) equation can be decomposed into the first two members 

of the coupled Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy by the binary nonlinearization of Lax pairs. 

In this paper, we construct the AT-th iterated Darboux transformation (DT) for the second- and third-order 

m-coupled AKNS systems. By using together the AT-th iterated DT and Cramer's rule, we find that the KPII 

equation has the unreduced multi-component Wronskian solution and the KPI equation admits a reduced 

multi-component Wronskian solution. In particular, based on the unreduced and reduced two-component 

Wronskians, we obtain two families of fully-resonant line-soliton solutions which contain arbitrary numbers 

of asymptotic solitons as y — ► =poo to the KPII equation, and the ordinary AT-soliton solution to the KPI 

equation. In addition, we find that the KPI line solitons propagating in parallel can exhibit the bound state 

at the moment of collision. 



1 Introduction 



In 1970, Kadomtsev and Petviashvili [1] derived the following (2+l)-dimensional nonlinear dispersive wave 
equation 



(1) 



{-Au t + 6uu x + u xxx ) x + 3a 2 u yy = 



±1), 



to study the stability of soliton solutions of the Korteweg-de Vries (KdV) equation with respect to weak 
transverse perturbations. Eq. (1) is usually called the Kadomtsev-Petviashvili (KP) equation, where a 2 = 
— 1 and 1 represent the KPI and KPII equations, respectively. As a natural extension of the KdV equa- 
tion in two dimensions, the KP equation (1) with both cases of a 2 has arisen in various physical contexts, 
e.g., plasma physics [1], surface and internal water waves [2], etc. The KP equation (1) is a universal 
completely-integrable (2+l)-dimensional nonlinear evolution equation (NLEE) [3] and its many integrable 
properties have been revealed in the past forty years, including the Lax representation [4], lump solutions 
to the KPI equation [5], line-soliton solutions to the KPI and KPII equations [6], an infinite number of 
symmetries and conservation laws [7], Painleve property [8], Backlund transformation [9], Darboux trans- 
formation (DT) [10], a recursion operator and bi-Hamiltonian structure [11, 12], and similarity reductions 
to the Painleve-type ordinary differential equations [14] and integrable (l+l)-dimensional NLEEs [15]. 

In recent years, the KPII equation has attracted intensive attention [16-23] because a large variety of 
soliton solutions of this equation were overlook in the literature before 2003. It is known that the KPII equa- 



tion admits the solution u 
form [24] 



(2) 



2tv-j In t(x, y, t) with the tau function t(x, y, t) expressible in the Wronskian 



T(x,y,t) 



h 

dx 



d N ~ 1 f 1 



h 

dh 
dx 



d N -\f 2 



dx r 



dx 1 



/JV 

df N 
dx 



dx N ^ 



where {fh\h=i are a set °f linearly-independent solutions of the linear system: af y = f xx (a = ±1), 
ft = fxxx- Ref. [17] has revealed a general family of multi-soliton solutions to the KPII equation by noting 
that the functions {fh\h=\ can be chosen as f h (x, y, t) = YaLi a h ie KlX+r7K ^ y+K f t (1 < h < N), where M 
is an positive integer greater than N, the phase parameters {ni}f£i are distinct, the constant coefficients a^i 
(1 < h < N; 1 < I < M) define an N x M coefficient matrix A := (a^), which is required to be of full 
rank [i.e., rank(A) = N] and all of whose nonzero N x N minors must be sign definite. 

If the matrix A is in the reduced row-echelon form (RREF) and satisfies the irreducibility conditions that 
each column of A contains at least one nonzero element and each row of A contains at least one nonzero 
element in addition to the pivot (first nonzero) entry, the function r in Eq. (2) can yield a multi-soliton 
configuration with N + = N asymptotic solitons as y — > +oo and AL = M — N asymptotic solitons as y — > 



— oo [17]. Moreover, it has been indicated [18-20, 22, 25] that the classification problem of soliton solutions 
to the KPII equation can be solved from the viewpoint of algebraic geometry and combinatorics because each 
equivalence class of (N, M — N) -soliton solutions corresponds to a derangement (which is a permutation 
with no fixed point) on {1, 2 . . . , M} with N excedances [19, 20], and each RREF N x M matrix A (all 
of whose maximal minors are non-negative) belongs to a unique totally non-negative (TNN) Grassmann 
cell which is a further decomposed cell of the Schubert cell of the TNN Grassmannian Gr^jf [25, 26]. In 
addition, the initial value problem of the KPII equation with V- and X-shape initial waves consisting of 
two distinct line solitons has also been studied by the direct numeric simulation [20-23], which gives an 
explanation of the Mach reflection phenomena in shallow water [20, 22]. 

Unlike the work in Refs. [16-23], the present paper is going to construct the multi-component Wronskian 
solution to the KP equation (1) and explore the underlying soliton solutions. The motivation of this study 
comes from the observation in two aspects. First, the binary nonlinearization of two symmetric Lax pairs 
can decompose the KP equation (1) into the first two members of the coupled Ablowitz-Kaup-Newell-Segur 
(AKNS) hierarchy [27, 28]: 



(3a) pj jV = -a 



_1 ( Pj,xx ~ 2 ^2 p n QnPj ) (1 < j < m), 
^ n=l ' 



(3b) qj sV = a 1 ( q j)XX -2^q n p n qj j (1 < j < m), 

^ n=l ' 

and 

m m 

(4a) pj :t = Pj.xxx ~ 3 ^2 P n ?™ Pj> ~ 3 ^2 Pn ' x qnP i ( X - J - m ^ ' 



ri=l n=l 



(4b) q jtt = qj,xxx -^^QnPn Qj,x ~ 3 ^ q n , x PnQj =0 (1 < j < 



_ m) 

n=l n=l 



where m is an arbitrary positive integer, Systems (3) and (4) are respectively referred to as the second- and 
third-order m-coupled AKNS systems, and 

m 

(5) u = -2^2pjqj, 

i=i 

exactly satisfies the KP equation (1). Second, the soliton solutions to both the vector nonlinear Schrodinger 
(NLS) equation and vector complex modified KdV equation, which are respectively two reductions of Sys- 
tems (3) and (4), have been found to be expressible in terms of the multi-component Wronskian [29-31]. 
Therefore, it is reasonable to infer that the KP equation (1) also admits the multi-component Wronskian 
solution. The structure of this paper is organized as follows. 



In Section 2, we follow the way in Ref. [32] to explicitly construct the iV-th iterated DT for Systems (3) 
and (4), which provides a basis for us to further obtain the multi-component Wronskian solution to the KP 



equation (1). We note that Ref. [32] has presented a general scheme of constructing the iV-th iterated DT 
for the AKNS hierarchy. However, there has been an absence of rigorous proof for such iV-th iterated DT. 
In this section, we give a complete proof that the Lax pairs associated with Systems (3) and (4) keep the 
form-invariance under the iV-th iterated DT. 

In most of the previous literature (see, for example, Refs. [24, 33]), researchers usually obtain the Wron- 
skian solutions by assuming that a given NLEE admits the Wronskian solution with certain condition, and 
then proving their conjecture by direct substitution the solution into the original equation or bilinear equa- 
tion^). In fact, one can also realize the Wronskian solutions to those Lax-integrable NLEEs by using 
together the iV-th iterated DT algorithm and Cramer's rule rather than relying on the guesswork [31]. More- 
over, if the Lax pair has been proved to be form-invariant under the DT, there is no need to verify the 
Wronskian solution again by means of the Pliicker relation. In Section 3, we derive the multi-component 
Wronskian solution to the KP equation (1) by combining the iV-fh iterated DT and Cramer's rule. 

In Section 4, we study the soliton solutions to the KP equation generated by the multi-component Wron- 
skian. In Section 4.1, we from an unreduced two-component Wronskian derive two families of fully- 
resonant line-soliton solutions to the KPII equation, which are shown to be two special cases of the soliton 
solutions generated by the Wronskian (2). In Section 4.2, based on a reduced two-component Wronskian, 
we obtain the ordinary iV-soliton solution to the KPI equation which, in general, describes the elastic col- 
lisions of line solitons. In particular, we find the bound state exhibited by the line solitons propagating in 
parallel at the moment of collision. 

Finally, in Section 5, we address the conclusions of this paper. 



2 N-th iterated Darboux transformation 



In the scheme of the (m + 1) x (m + 1) -matrix AKNS scattering formulation [34], the Lax representation 
of Systems (3) and (4) can be written in the form 

(6a) $ x = U(X)<f> = (\U + U 1 )<f>, 

(6b) % = F(A)$ = (X 2 V + A Vi + V 2 ) 

(6c) <$> t = W{\)<& = (A 3 Wo + >?Wi + XW 2 + W 3 ) <£, 

with 



-E m 




Vo = --U , Vt = 

a a a I „t _r,T 



q p, 



W = 4 U , W! = 4 Ui, W 2 = -2aV 2 



p xx - 2pq T p 



pqj - Pxq T 



ilx - 2 q T pq T q T Px - qjp 



where A is the spectral parameter, = (cf>i, (f>2, • • • , 4>m.+i) T C denotes the transpose of vector) is the 
vector eigenfunction, p = (pi,P2, ■ ■ ■ ,Pm) and q = (q\, q 2 , . . . , q m ) are the vector potentials, E m is the 
m x m identity matrix, and the compatibility conditions $> xy = $> yx and & xt = <& t x yield Systems (3) 
and (4), respectively. 

According to the idea of the DT [32], it requires that under the iV-th iterated eigenfunction transformation 

(7) $' = T(A)$, 



,P'm) 



the new eigenfunction & = <fi 2 , ■ ■ ■ , 4>' m +i) a l so satisfies System (6) with the matrices U(X), V(X) 
and W(X) replaced respectively by U'(X), V'(X) and W'(X) in which the new potentials p' = {p'i,p' 2 
and q' = (q[, q' 2 , . . . , q' m ) are respectively instead of the old potentials p and q, that is, 

(8) & x = U'(\)&, & y = V'(X)&, & t = W'(\)&. 

Thus, the Darboux matrix T(A) has to obey the following three conditions: 

(9a) T X (X) + T(X)U(X) = U'(X)T(X), 

(%) T y (X) + T(X)V(X) = V'(X)T(X), 

(9c) T t (A) + T(X)W(X) = W'(X)T(X), 

where Conditions (9a)-(9c) respectively coiTespond to the form-invariance of Eqs. (6a)-(6c). 

For the satisfaction of Conditions (9a)-(9c), we assume that the matrix T(A) is of the form [32] 

( a(X) A (A) ••• (3 m (X)\ 



(10) 

with 
(11) 

(12) 

(13) 



T(A) 



7i(A) S n (X) 

\lm(X) 5 ml (X) 



Slm(X) 
u rnrn 



N-l 



N-l 



a(X) = X N -J>Wa», Pj(X) = £ ^\-X) n (1 < j < 



m) 



n=0 



n=0 



N-l 



7i(A) 



^ 7i (n) A" (l<t<m), fe(A) = A 7V +^4 n) (-Af (l<i<m), 



JV-l 



(n). 



n=0 
N-l 



n=0 



Sij(X) = 4f(-^) n (1 < i,3 < m; 



n=0 



where (3- ,7^ and 8\a (1 < i, j < m; < n < N — 1) are the functions of y and i which are 
determined by 

(14) T(X k )$ k = (1 < A; < K = mN + N), 

with = (fk,9u , ■ ■ ■ jS'fc ) T as the solution of System (6) corresponding to A = A& (A& 7^ Af, 1 < 
k, I < K). It should be noted that {^ k } k=1 we a set °f linearly-independent solutions of System (6) because 
Afc 7^ A^. Accordingly, one can uniquely determine P^, 7^ and 5^ (1 < i,j < m; < n < N — 1) 
from (14) which can be expanded as 

m m 

(15) a(A fc )+^^' ) ^(A fe ) = 0, li (X k ) + J2^k )s ij( X k) = (1 < i < m; 1 < k < K), 

j=i 3=1 

where w[ j) = g[ j) /f k (1 < j < m; 1 < k < K). 

Before proving that Conditions (9a)-(9c) are satisfied with the Darboux matrix T(A) given by Eq. (10), 
in which /3j n \ 7}^ and 8^ (1 < i, j < m; < n < N — 1) are determined by Eqs. (15), we present the 
following two important lemmas. 

Lemma 2.1. Let the Darboux matrix T(X) be in the form of (10) with /3j , 7^ and 8^ (1 < i,j < m; 
< n < N — 1) determined by Eqs. (15). Then, the determinant ofT(X) can be expanded as 

K 

(16) detT(A) = [](A-A fc ), 

fc=i 

where X k (1 < k < K) are the roots o/det T(A). 

Proof. It is easy to know that detT(A) is a monic polynomial of degree K. On the other hand, one can 
employ Eqs. (15) to express the first column of T(X k ) as the linear combination of other columns: 

m 

(17) [a(A fc ),7i(A fc ), • • • ,7m(A fc )] T = - J^^lPjiWij (A fe ), . . . , 5 mj (X k )] T (l<k<K), 

3=1 

which suggests that X k (1 < k < K) are the roots of detT(A). Therefore, the determinant of T(A) can be 
expressed in the form of (16). □ 

Lemma 2.2. Let us define that 

(18a) [u h i(X)] {m+l)x{m+1) = [T X (X) + T(A)*7(A)]T*(A), 

(18b) [^(A)] (m+ i)x (m +i) = [T y (X) + T(X)V(X)]T*(X), 

(18c) [%(A)] (m+ i)x(m+i) = [T t (X)+T(X)W(X)]T*(X), 



where T*(A) is the adjoint matrix ofT(X). Then, X k (1 < k < K) are the roots of Uhi(X), Vhi(X) and 
w hi(^) (1 < h,l < m + 1), that is, 

(19) u hl (X k ) = 0, v M (X k ) = 0, w hl (X k ) = (1 < h, I < m + 1; 1 < k < K). 



The proof of Lemma 2.2 is given in Appendix A. In the following, we use deg[/(A)] to represent the 
degree of the polynomial /(A), and a* (A), /3*(A), 7* (A) and 5*j(X) to denote the algebraic cofactors of 
a(A), /3j(A), 7i(A) and Sij(X) (1 <i,j< m), respectively. 

Proposition 2.3. Assume that Q k = (fk,9^\ ■ ■ ■ > dk"^)^ satisfies Eq. (6a) with X = X k , where 1 < k < K . 
Then, the Darboux matrix T(A) given by (10) obeys Condition (9a), provided that pp , 7^ and 5^ 
(1 < i, j < m; < n < N — l)are determined by Eqs. (15), and the N-th iterated potential transformations 
are given by 

(20) p / = p-2(-l) JV - 1 b( JV - 1 ), 

where = and = (jf^, . . . 



Proof. Let 11(A) = [T X (X) + T(A)J7(A)]T*(A)/detT(A). The calculation of algebraic cofactors a*(X), 
/3*(A), 7 *(A) and tff.(A) yields that deg[a*(A)] = deg[^(A)] = mN and deg[/3*(A)] = deg[ 7 *(A)] = 
deg[<5*j(A)] = mN—1 (1 < i, j <m;i^ j), which tells us that deg[ti/ lfe ( A)] = K+l and deg[u w (A)] = K 
(1 < h, I < m + 1; h ^ I). On the other hand, Lemmas 2.1 and 2.2 imply that u/^(A) (1 < h, I < m + 1) 
can be exactly divided by det T(A). Therefore, the matrix n(A) can be written in the form 

\ 



Ml 



(21) 



n(A) 



7T 



+ 4? 

(0) 



7T 



(0) 
12 



7T 



21 



A 7T22 + vr. 



(0) 
22 



(0) 

l,m+l 

7T (0) 
^m+l 



7T 



(0) 

m+1,1 



vr (0) 
^+1,2 



X (!) 1 (0) 

A ^m+l.m+l ' ^m+l.m+l / 



where it}) and 7r^ ; (1 < /i, I < m + 1) are all the functions dependent on x, y and t. By expanding 
T X (X) + T(A)f7(A) = n(A)T(A) and comparing the coefficients of X N+l and A^, we can find 



7T 



7T 



(1) _ 
11 — 

(0) 



vr (0) 
j ^11 



0, 

-1) 



7T (1) - 

N-la(N-l) 
Pj ' 



7T 



(0) 



7T 



(0) 

t+1,1 



(l<i,j<m), 
27f 7V " 1) (l<i,j<m). 



It is obvious that the matrix n(A) is exactly equal to U'(X) with the new potentials p' and q' given by (20), 
which means that T X (X) + T(X)U(X) = U'(X)T(X). □ 



Similarly, we present other two propositions, as stated in the following: 



Proposition 2.4. Assume that &k = gjp, . . . , g^) T satisfies Eq. (6b) with X = where 1 < k < K. 
Then, the Darboux matrix T(X) in (10) obeys Condition (9b), provided that (3j n \ 7^ and 5^' (1 < i,j < 
m; < n < N — 1) are determined by Eqs. (15), and the new potentials p' and q' are given by (20). 

Proposition 2.5. Assume that = (/^ , > • • • > 9k^ ) T satisfies Eq. (6c) with X = X^, where 1 < k < K. 
Then, the Darboux matrix T(X) in (10) obeys Condition (9c), provided that (3j, 7}"^ and 5^ (1 <i,j < 
m; < n < N — 1) are determined by Eqs. (15), and the new potentials p' and q' are given by (20). 

Seen from Propositions 2.3-2.5, the Darboux matrix T(X) makes sure that the new eigenfunction <£' = 
T(A) $ also satisfies Lax representation (6) for the new potentials p' and q' given by (20). That is to say, 
the compatibility conditions & xy = & yx and & xt = & tx give the same systems (3) and (4) except for p' and 
q' instead of p and q, respectively. Therefore, we come to the following theorem: 

Theorem 2.6. Suppose that {&k}k=i are a set of linearly-independent solutions of System (6) which corre- 
spond to a set of different spectral parameters {Xk}^ =1 - The eigenfunction transformation (7) and the po- 
tential transformations (20) constitute the N-th iterated DT (<£, q, p) — > ($', q', p') of Systems (3) and (4), 
where the Darboux matrix T(X) is determined by Eqs. (15). 

Remark 2.7. With some constraints between the potentials p and q, Systems (3) and (4) can be reduced 
to the NLEEs belonging to some known integrable hierarchies such as the coupled NLS hierarchy [35 ] 
(p = ±q, where the bar represents complex conjugate) and coupled modified KdV hierarchy [36] (p = ±qj. 
In the above way, one can also use Transformations (7) and (20) to constitute the N-th iterated DTfor those 
reduced cases, but {<&k\k=i are reau ^ re d to be particularly selected. For example, ifp = ±q, the functions 

{^k\k=l are cnosen as [31] 

(22a) $( fe -i)( m+ i) + i = (/ fc ,5i 1) ,...,5i m) ) T (1<*<JV), 

(22b) * (fc _ 1)(m+1)+i+1 = (±^ ) ,0~^0,/ fc ,0^~0) T (1 < k < N; 1 <j < m), 

where ^^-i)(m+i)+i satisfies System (6) with p = ±q and X = (1 < k < N), and ^^-i)( m +i)+j+i 
( 1 < j < m) are orthogonal to $(fc_!)( m+1 ) +1 . 

3 Multi-component Wronskian solution 

In this section, by combining the iV-th iterated DT of Systems (3) and (4) and integrable decomposition from 
the KP equation (1) to Systems (3) and (4), we try to construct the multi-component Wronskian solution to 
the KP equation (1). 



In order to solve the functions $ N ^ and 7^ ^ (1 < i,j < m), we rewrite Eqs. (15) in the following 



form 
(23) 
(24) 
where 



-G 



(m) 



X=(Af/i, 



{F K 



r (i) 



Fit 



xN 



( fi A1/1 •• 
\ /a A#/a- • ■ 



G 



0') 

AxA 



(J) 



V <?a "A as 



A' 



/nr(m) \ v _ (\N(i) 
~ Lj KxN) Y i - l A l 9l ; 

Af-Vi \ 

Arv 2 



A A / K J 

■ (-A 2 )^ 1 5 ? ) 
•• ("Aa)^ 1 ^, 



A^/a) T , 

■XkSV (1<<< 



m 



(1 < j < m), 



X — (a^ \ . . . , ^ ; f3^ \ f3^ \ . . . , f3^ ; . . . ; /3^\ . . . , ^) 



(°) 



7i .7* 



l )•••)/! 



(iV-1) ,(0) ,(1) 



' °il ' °il ' ■ ■ ■ ' °il ' • • • ' °im > °im ' • • • > °: 



r (tf-l)\T 



) (1 < i < m). 



By employing Cramer's rule, we can obtain the functions pf 1 ^ and 7^ ^ as 



(25) 

with 

(26a) 
(26b) 
(26c) 



p (N-l) = ( _ 1)jN -t^ 7 (iV-l) =( _ 1) (^l)A-lX i _ (1 < iti < ro)j 



(1) 



(2) 



t = det(F K xN, -G 



(i) 

AxiV> ' ' ' ' 



-G 



(m) 
KxN 



*3 



(1) 



V, 



(2) 



det(Fj 
det(F j 



Ax(7V+l) 



,-G 



Ax(A-l), KxNi 



G 



(1) 

AxW' 
(1) 



,-G 
,-G 



Cj) 



G 



(m) 



Ax(A-l)' ' "-'KxN 



) (1<J< 



m 



(0 

Ax(A+l)' 



,-Gg^) (l<t<m). 



where r, xf 1 andx^ (1 < i,j < m) are the multi-component Wronskians which have been proposed in 
Ref. [29]. Thus, the iV-th iterated potential transformations (20) can be written as 



(27) 



+1)N+1 

T 



(1) 



.(2) 



2 (_i)(i-i)*-iXL_ (l<i,j<m). 



With p = q = 0, the general solution of System (6) with A = A& (A& 7^ A;; 1 < k, I < K) is given as 
follows: 



(28) 



= {fk,9k\---,9k 



(a fc e*S&<V***, . . . ,6j™e-3"*) 1 (1 < fc < tf), 



where 
(29) 



e k = 2(\ k x-2a\ 2 k y + 4\lt). 



It is noted that the field u of Eq. (1) is usually assumed to be a real one. For the case a 2 = 1 which 
corresponds to the KPII equation, we require the parameters X k , a k and &i (1 < k < K; 1 < j < m) be 
all real constants. 

Proposition 3.1. With Q k = (f k , g£\ . . . , given by (28) for 1 < k < K, the KPII equation 

admits the solution u = 2 (lnr) xa; , where the tau function r is expressed as the following unreduced multi- 
component Wronskian: 



(30) 
with 



T — \AkxK@KxK^-KxNi ^KxK^KxK^KxNi • • • > 



'KxK KJ KxK lv KxN\i 



A K xK = diagfai, 



,a K ), B ij) 



KxK 



diag( 



di) 

h > ■ ■ ■ j 



(1<J< 



Q+ xK = diag(e5 e \ . . . ,e5 9 «), Q^xK = diag(e~^ 



A + 

KxN 



( 1 Ai 
1 A 2 

^1 A A - 



N-l\ 



,7V-1 



A 



(l 
1 



7V-1 



-A-i 
-A 2 



(-A 2 ) JV - 1 



Proof. With p = q = 0, substituting Eqs. (27) into Eq. (5) yields 
(31) u- 



m m v ( 1 ) v ( 2 ) 

2E^ = 8 E^ 



On the other hand, one can follow part (ii) of Theorem 3.3 in Ref. [29] to obtain the following multi- 
component Wronskian identity: 

1 , 



(32) 



y- (i) (2) 



1" T » ~ T x) > 



which is substituted into the right-hand side of Eq. (31), giving that u = 2 (In r) xx . □ 



In the case a 2 = —1, the solution u = 2 (lnr) xx with r given by (30) also satisfies the KPI equation, 
but it does not meet the requirement that u must be a real field. According to Remark 2.7, if imposing the 
constraint p = eq (e = ±1) on Systems (3) and (4), one can choose {& k }f =1 as 

(33a) $ (fc _ 1)(m+1)+1 = (a fc e^4VH,...,4 m ) e -^) T (1 < k < N), 

j-l m-j 

(33b) $ {fc _ 1)(m+1)+j+1 = (4 i) e-^"^07^T0,aA : e^"^0~^T0) T (1 < k < N; 1 < j < m), 
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where the phase 9k is the same as given in Eq. (29), but the parameters and or (1 < k < N; 

1 < j < m) are all complex constants. Then, p'- and its complex conjugate p 1 - are obtained as follows: 
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where the functions r, Xj an d Xj we given by 
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where F NxM = A NxN & NxN „ NxM , - v . A/ 
multi-component Wronskian identity (32), we know that 



B NxN®NxN A NxM (M = N-1,N,N + 1). Via the 
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which implies that u is a real function. 



5fc"^) T gi yen by Eqs. (33a) and (33b) for 1 < k < K, the 



Proposition 3.2. With = (fk,g k 
KPI equation admits the solution u = 2 (lnr) xx , where the tau function r is a reduced multi-component 
Wronskian given in Eq. (35). 



4 Soliton solutions to the KP equation 



In this section, we explore the soliton solutions to the KPII equation generated by the unreduced multi- 
component Wronskian (30), and to the KPI equation generated by the reduced multi-component Wron- 
skian (35). 



4.1 Fully-resonant soliton solutions to the KPII equation 



As for the KPII equation, one might naturally ask whether the soliton solutions generated by the multi- 
component Wronskian (30) are different from those by the single Wronskian (2). Since the multi-component 
Wronskian (30) with m > 1 generates singular solutions for the generic choice of parameters, we in this 
subsection discuss the soliton solutions in u = 2 (In r) xx with 

(39) T = \A2Nx2N<3>2Nx2N^2NxN> ^2Nx2N®2Nx2N^2Nxn\- 

For convenience, we take the coefficient matrix B^ x2N = B2NX2N = diag(— &i, . . . , — &2iv) with b k = 
for k E [2iV] := {1,2,..., 2iV}, and for comparison with the soliton solutions in Refs. [17-20] we 
assume that 9^ = —K k x — crn^y — k\H, where K k = — 2\ k and {Kk} k =i are without loss of generality 
ordered as k\ < k,2 ■ ■ ■ < K2N- 

By using the Laplace expansion technique and Binet-Cauchy theorem, we can expand the function r in 
Eq. (39) as follows: 
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where In = {h, ■ ■ ■ ,i/v} an d Jn = {ji, ■ ■ ■ ,Jn} are two subsets of [2N] with 1 < i\ < ■ ■ ■ < ijy < 2N 
and 1 < ji < • • • < j'jv < 2N. In order to ensure that the solution u to the KPII equation resulting from the 
function r in Eq. (39) is non-trivial and non-singular, the coefficients a±, . . . , a,2N, &i, • • • , ^2Af are required 
to satisfy the following three conditions: 

(i) afc and bk are not equal to zero for the same k G [2N] ; 

(ii) \{a k \a k / 0}| = N + L and \{b k \b k / 0}| = iV + M, where 1 < L, M < N and | • | denotes the 
number of elements in a set; 

(iii) a k a k+ ib k b k+ \ < for 1 < k < 2A^ — 1. 



Here, condition (i) requires that the function r does not reduce to zero; condition (ii) requires that the 
function r contains at least two exponentials; condition (iii) requires that the function r is sign definite, so 
that r has no zeros for all (x, y, t) € M 3 . 



Considering that there might be some zeros among the coefficients a±, . . . , a,2N, &i> • • • , ^2N, we assume 
that Of. = for n € [£'] and fy. = for n 6 [M'], where L' = iV - L, M' = iV - M, 1 < z* < ■ ■ ■ < 
i* L , < 2N, 1 < j* < ■ ■ ■ < j* M , < 2N and i* / j*. Thus, the function r in Eq. (40) can be expressed as 
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where Z A / = {n, . . .,i M }, Jl = {ji, ■ ■ -Jl}, X* v = {**, • • ■ Jt v = {jl, ■ ■ ■ ,j* M ,}, T** and Tj* 
are given by 
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With scaling transformation, the function r in Eq. (41) is further equivalent to 
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in generating the solution u of the KPII equation. Without loss of generality, we suppose that Xm U = 
[M + L] and X* L , U J^, = {M + L + 1, . . . , 2N}. We notice that the function r' in Eq. (43) contains 
all possible combinations of M phases out of {9\, . . . , 9m+l}> which is the same as the case of a family 
of fully-resonant soliton solutions to the KPII equation discussed in Ref. [16]. Therefore, we immediately 
come to the following results: 

Proposition 4.1. The dominant exponentials of the function t 1 in adjacent regions of the xy-plane as ay —¥ 
±00 contain M — 1 common phases and differ by only one phase, and the transition between two such 



exponentials occurs along the line defined by Ly : B{ = 9j where j = i + M for i G [L] as ay — > oo and 
j = i + L for i 6 [M] as ay —> — oo. In a neighborhood of the transition line Ly as ay — > ±00, the 
asymptotic behavior of u = 2 (In 
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where defines the asymptotic line soliton as ay — > ±00, and A^jj, ^-[ij] and Q[ij] respectively 
represent the amplitude, wave vector and frequency of asymptotic soliton. 

We remark that the function t' in Eq. (43) generates two families of fully-resonant soliton solutions: 
one is (L, M)-soliton configuration with L asymptotic solitons as y — >■ 00 and M asymptotic solitons as 
y — > —00 for a = 1 (see Figure 1), and the other is (M, L)-soliton configuration with M asymptotic solitons 
as y — > 00 and L asymptotic solitons as y — > —00 for a = — 1 (see Figure 2). Such two families of solutions 
correspond to two special cases of the soliton solutions generated by the function r in Eq. (2) because the 
soliton amplitude Aha, wave vector Kuj-i and frequency Qua in the asymptotic expression (44) are the 
same as those given in Refs. [17, 19, 20]. 



4.2 Ordinary iV-soliton solution to the KPI equation 



In this subsection, we are devoted to discussing the soliton solutions to the KPI equation generated by the 
function r in Eq. (35) with m = 1 and e = — 1, i.e., 
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where r is a real- valued function and has no zeros for all (x,y,t) € M 3 . Our consideration is based on 
the following two facts: (a) The solutions generated by the (m + l)-component Wronskian (35) do not 
contain more valid free parameters than those by the two-component Wronskian (45); (b) The function r 
in Eq. (35) with e = 1 has zeros for some (x, y, t) G IR 3 . For convenience, we define that = and 
= \ (fJ-k + f^fc) for k € [N], where a = ±i, fik and V}. are real constants. Without loss of generality, we 
assume that a\. = 1, / for k € [N] and are well ordered as v\ < < ■ ■ ■ < vn- 
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Figure 1. A fully-resonant (3, 2)- 
soliton solution to the KPII equa- 
tion generated by the function r' in 
Eq. (43) with the parameters cho- 
sen as L = 3, M = 2, a = 1, 



h' — h' — h' — h 1 — h' - 
°1 — °2 — °3 — °4 — °5 — 

Kl = —0.8, K2 = —0.35, K3 

0.25, k 4 = 0.65 and k 5 = 1.25. 



1. 



Figure 2. A fully-resonant (2,3)- 
soliton solution to the KPII equa- 
tion generated by the function r' 
in Eq. (43) with the same parame- 
ters as those in Figure 1 except for 
a = -1. 



For the case N 
(46) 



1, the solution u = 2(mr) 
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u 



/^sech 2 



with r given by Eq. (45) can be expressed as 

In |6x| 



1 1 1 / 2 

L2 Mia; ~ 2^ lVlV+ S^ 1 ^ 1 



3z/i) t 



which describes the one-soliton solution to the KPI equation (see Figure 3), where the amplitude, wave 
vector and frequency of soliton are respectively given by | [if, fii, — | ii\v{) and 1 1^1(1^1 — 3isf). It is 
mentioned that such line-soliton solution is in accordance with those obtained by the Hirota method [6] and 
inverse scattering transform [37]. 



For the case N > 2, we can employ the results of Ref. [29] to analyze the asymptotic behavior of the 
function r in Eq. (45). According to Lemma 4. 1 in Ref. [29], we know that the function r in Eq. (45) can be 



expanded as the sum of such exponential terms as exp [| J2k=i { r kGk + Sfc#fc) ] . where r^, Sk G {—1,1}, 
|{^fc|^fc = — 1}| = |{s/c|sfc = — 1}| and |{rfc|rfc = 1}| = = 1}|. Further using Theorem 4.4 in 

Ref. [29], we obtain that along the line 6 n + 9 n = (1 < n < N) as y — > ±oo for finite values of t, the 
solution u = 2{\ut) xx is asymptotically determined by 



u ~ u, -, 



+ r,t ^P { ^ [ - 9 n - n + £4(0* + 9k)] }) 



= ^2 sech 2^^2_Ji +ln ^ c ± /?? ± J (l< n <JV), 



where Ur n defines the n-fh asymptotic line soli ton [n,n] as y — > ±oo, and 77* correspond to the 
coefficients of two dominant exponentials in the expansion of r, the parameters e^ n 's are defined as 

__/-!, for k G UB [ n\ + _ J 1, for fc G B^ U b1 E) , 
^ " \ 1, for A; G B^ U Bf\ " \ -1, for k G B« U B^, 

withfij? = > 0,/ = l,...,n- 1}, b£° = {l\tn < 0,1 = n + l,...,N}, B n W) = {l\m < 0,1 = 
l,...,n-l}andB^ = {l\m > 0, 1 = n + 1, . . . , N}. 




Figure 3. A one-soliton solution 
to the KPI equation generated by 
the function r in Eq. (45) with the 
parameters chosen as N = 1, a = 
i, fi\ = 1, vi = 2 and bi = 1. 



Figure 4. An obliquely-colliding 
two-soliton solution to the KPI 
equation generated by the function 
r in Eq. (45) with the parameters 
chosen as N = 2, a = i, \i\ = 
H2 = 1, v\ = —1.6, V2 = 2 and 
h = b 2 = 1. 



The asymptotic expression (47) shows that the n-th asymptotic soliton as y — > oo has the same velocity 
and amplitude as the n-th asymptotic soliton as y — > — oo except for a slight phase shift, as displayed in 
Figure 4. Therefore, the iV-soliton solution generated by the function r in Eq. (45) reflect only the ordinary 
elastic collisions of line solitons in the KPI equation. Particularly with v n = u^{n ^ k), the n-th and fc-th 
asymptotic solitons propagate parallel to each other (see Figure 5), but they exhibit the bound state [38] at 
the moment of collision (see Figure 6). It is an interesting question whether the KPI equation also admits 
the inelastic soliton collisions (in which the numbers, amplitudes or directions of asymptotic solitons are not 
the same as y — > ±oo), like the cases in the KPII equation. 
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Figure 5. Parallel propagation of 
two line solitons generated by the 
function r in Eq. (45) with the pa- 
rameters chosen as N = 2, a = i, 
Hi = 1, u 2 = 1.5, v\ = U2 = 0.8 
and 6i = 62 = 1. 
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Figure 6. Bound state between 
two line solitons at the moment of 
collision generated by the function 
r in Eq. (45) with the same param- 
eters as those in Figure 5. 



5 Conclusions 



In this paper, we have constructed the iV-th iterated DT for the second- and third-order m-coupled AKNS 
systems (3) and (4), and rigorously proved the form-invariance of their Lax representation (6). Then, using 
Cramer's rule and based on the integrable decomposition from the KP equation (1) to Systems (3) and (4), 
we have derived the unreduced multi-component Wronskian solution to the KPII equation and a reduced 
multi-component Wronskian solution to the KPI equation, which implies that the multi-component Wron- 
skian is a third form of the tau function of the KP equation in addition to the Wronskian and Grammian. 
Further, by analyzing algebraic properties of the unreduced and reduced two-component Wronskians, we 



have obtained two families of fully-resonant line-soliton solutions which, in general, contain unequal num- 
bers of asymptotic solitons as y — > =foo to the KPII equation, and the ordinary iV-soliton solution which 
generally describes the elastic collisions of line solitons to the KPI equation. In particular, we have found 
that the KPI line solitons which propagate in parallel can exhibit the bound state at the moment of collision. 
It is remarkable that the double Wronskian representation of the iV-soliton solution to the KPI equation has 
been revealed in this work for the first time. 
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Appendix A Proof of Lemma 2.2 



Proof. First, we make use of Eqs. (15) to compute T*(X k ) as 
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(1 < i < m). 



Then, we expand [T^A*) + T(A fc )C/(A fc )]T*(A fe ), [r y (A fc ) + T(A fc )y(A fc )]T*(A fc ) and [T t (X k ) + T(A fc ) 
VF(A fc )]T*(A fc ), obtaining that 
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Using Eqs. (15) and recalling that = (f k ,g9 , . . . ,g["^) T satisfies System (6) with A = X k , we 
obtain the derivatives of a(X k ) and 7i(Afe) with respect to x, y and t as follows: 
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By means of Eqs. (15) and (A.8)-(A.13), we remove ct(A fe ), 7i(A fc ), a x (A fc ), 7i, x (A fc ), a y (A fc ), 7^ 
ctt(\k) and 7i,t(Afc) from the right-hand sides of Eqs. (A.2)-(A.7), and finally prove that Uf l i(X] c ) 
Vhii^k) = and w h i{X k ) = (1 < h, I < m + 1). □ 
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